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Abstract 

The study of "Entanglement of Formation" of a mixed state of a bipartite system in 
high-dimensional Hilbert space is not easy in general. So, we focus on determining the 
amount of entanglement for a bipartite mixed state based on the concept of decomposable 
optimal entanglement witness (DOEW), that can be calculated as a minimum distance 
of an entangled state from the edge of positive partial transpose (PPT) states which has 
the most negative (positive) expectation value for non-PPT (bound) entangled states. 
We have constructed DOEWs based on the convex optimization method, then by using 
of it we quantify the degree of entanglement for two spin half particles under the Lorentz 
transformations. For convenience, we restrict ourselves to 2D momentum subspace and 
under this constraint when the momentum and the Lorentz boost are parallel, we have 
shown that the entanglement is not relativistic invariant. 

Keywords : Relativistic entanglement, Measure of entanglement, Optimal 
entanglement witnesses, Convex optimization 
PACS numbers: 03.67.Hk, 03.65.Ta 

1 Introduction 

In the recent years it became clear that quantum entanglement p] is one of the most important 
resources in the rapidly growing field of quantum information processing because, quantum 
entangled states produce nonclassical phenomena. Therefore, specifying that a particular 
quantum state is entangled or not, must be important since for separable quantum states 
statistical properties can be explained entirely by classical statistics. 

If a density matrix, P(a,b), of a composite system (A,B) can be written as a sum of products of 
density matrices of its components, p(A) and P(b) in the form P(a,b) = * P\a) ® P\b)i ® ^ 
qi ^ 1 and J2i Qi — 1j then the system is separable, otherwise it is entangled. The first and 
most widely used related criterion for distinguishing entangled states from separable ones, is 
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PPT criterion, introduced by Peres [2]. Furthermore, a necessary and sufficient condition for 
separability in Ti 2 <8> H-2 and 7i 2 ® Hilbert spaces ( TCd denotes the Hilbert space with 
dimension d endowed with usual inner product denoted by ( . )) was shown by Horodecki 
P], which was based on a previous work by Woronowicz [3]. However, in higher dimensions, 
there are PPT states that are nonetheless entangled, as was first shown in Ref. [5], again 
based on Ref. [14]. These states are called bound entangled states because they have the 
peculiar property that no entanglement can be distilled from them by local operations [5]. So 
the PPT criterion is not sufficient for separability. Another approach to distinguish separable 
states from entangled ones is entanglement witness (EW) [11]. An EW for a given entangled 
state p is an observable W, whose expectation value over all separable states is nonnegative, 
but strictly negative on p. There is a correspondence between EWs and linear positive (but 
not completely positive) maps via Jamiolkowski isomorphism [12J. As an example the partial 
transposition is a positive map (PM). Despite of the fact that EWs are designed mainly for 
detection of entanglement, it has been shown [13] that the optimal EW associated with a 
density matrix p -in a the sense that, the expectation value of the optimal EW (associated 
with p) over p is the most negative value between the expectation values of other EWs over p- 
can be used as measure of entanglement quantifying the amount of entanglement of p. 

Relativistic aspects of quantum mechanics have recently attracted much attention, espe- 
cially in the context of the theory of quantum information and entanglement. Recently, several 
groups have focused their investigations on relativistic quantum entanglement [171 [18j [19j [20], 
HH [221 1231 [2H [251 12SI E7J EH [291 EH E21 S3]. Peres et al. [33] have observed that the reduced 
spin density matrix of a single spin-| particle is not a relativistic invariant, given that Wigner 
rotations entangle the spin with the particle momentum distribution when observed in a mov- 
ing frame [35]. Gingrich and Adami have shown that the entanglement between the spins of 
two particles is carried over to the entanglement between the momenta of the particles by the 
Wigner rotation, even though the entanglement of the entire system is Lorentz invariant [22] . 



4 



Terashimo and Ueda [23] and Czarchor [21] suggested that the degree of violation of the Bell 
inequality depends on the velocity of the pair of spin-| particles or the observer with respect 
to the laboratory. Alsing and Milburn studied the Lorentz invariance of entanglement, and 
showed that the entanglement fidelity of the bipartite state is preserved explicitly. Instead of 
state vector in Hilbert space, they have used a 4-component Dirac spinor or a polarization 
vector in favor of quantum field theory [25]. Ahn also calculated the degree of violation of 
the Bell's inequality which decreases with increasing of velocity of the observer [26]. Most 
of the previous works were concerned with pure states although authors in [271 1281 129] have 
considered mixed quantum states that are described by superposition of momenta. In pervious 
works the measure of entanglement for pure states under the Lorentz transformation has been 
calculated using the Bell's inequality. 

In Ref . |30j . we investigated spin-momentum correlation in single-particle and showed that 
entanglement decreases with respect to the increasing of observer's velocity. However, in this 
work, we are concerned with degree of entanglement for two relativistic spin half particles and 
will show that it is possible to use the idea of DOEW operators as measure of entanglement 
in high-dimensional Hilbert space because finding the region of separable states is not easy in 
general. Instead of using a superposition of momenta for particles we use only two momentum 
eigen states {jp\ and p?) and restrict ourselves to 2D momentum subspace and under this 
constraint, our analysis will concentrate to indicate a new method for constructing DOEWs 
based on the convex optimization method which has been widely used in quantum information 
theory [7J El O HH] • Likewise, in Refs. [361 [37] a connection between Hilbert-Schmidt measure 
and optimal EW associated with a state has been discussed. However, we show that the 
constructed DOEW is based on minimal distance of the corresponding entangled state from 
the edge of PPT states. So, the minimal distance or expectation value is most negative 
(positive) for non-PPT (bound) entangled states. At the end, by using the DOEW, we show 
that when momentum and Lorentz boost are parallel, the entanglement of the entire system 
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is not invariant under the Lorentz transformation. 

This paper is organized as follows: Sec. II, is devoted to two-particle relativistic quantum 
states. In Sec. Ill, we present the DOEWs which are constructed using the convex optimization 
method. In Sec. IV, we derive a simple formula for optimality of the decomposable EW via 
Hilbert- Schmidt measure. In Sec. V, we explicitly calculate the amount of entanglement of 
two relativistic particles. The last section contains concluding remarks. 

2 Two relativistic particles quantum states 

Suppose we have a bipartite system with its quantum degrees of freedom distributed among 
two parties A and B with Hilbert spaces Ha and Hb, respectively. The quantum states of 
one particle is made by two degrees of freedom: momentum p and spin. The former is a 
continuous variable with Hilbert space of infinite dimension but in this work, we have only two 
eigen states pi and p 2 , while the latter is a discrete one with Hilbert space of spin particle. 
The pure quantum state of such a system can always be written as 

2 -1/2 
i=l j=l/2 

where both \pi) and \p 2 ) are two momentum eigen states of each particle and the kets ||) 
and 1^) are the eigenvectors of spin operator a z . c^'s are complex coefficients such that 
Ylij \ c ij\ 2 = 1- Using the general pure state (12.11) . we define four orthogonal maximal entangled 
Bell states for one-particle in terms of momentum and spin states as follows, 

M 

|V>3> 



= ^(|Pl> «> |g> + |P2> «) | — -», 
= ^{\Pl) ®\\) ~\P2) ®\~\)), 
= ^{\P2)®\\) + \Pl)®\-\)), 
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|^4> = ^(b 2 >®|^>-bi>®l-^>)- 



(2.2) 



Then, using the upper Bell states, we construct two-particle quantum states as 



IV>±> (1 ' 2) = 


1 


IV>i)I^i>± 


1^2)1^2)), 


|0 ± ) (1 ' 2) = 


1 


1^1)1^2) ± 




l^) (3 ' 4) = 


1 


1^3)1^3) ± 


1^4)1^4)), 


|0 ± ) (3 ' 4) = 


1 


1^3)1^4) ± 


1^4)1^3)) 


l^) (1 ' 3) = 


1 




l^s) 1^3)), 


|0 ± ) (1 ' 3) = 


1 


|V>1>N>3>± 


l^3)|^l)) 


l^) (2 ' 4) = 


l 


1^2)1^2) ± 


1^4)1^4)), 


|0 ± ) (2 ' 4) = 


1 

7i ( 


1^2)1^4) ± 


1^4)1^2)) 


l^) (1 ' 4) = 


1 

7t 




1^4)1^4)), 


|0 ± > (M) = 


1 

7! ( 


l^l)l^4) ± 





l<M (2,3) = ± M*)l&», |0 ± ) (2 ' 3) = ± (2.3) 

and introduce the following orthonormal entangled states in 16- dimensional Hilbert space: 



|$!) = cos0|<0+> (1 ' 2) + sin^ + ) (3 < 4) , |$ 2 > = cos6>|^_> (1 ' 2) + sin^|^-) (3 ' 4) , 

|$ 3 > = cos^ + >( 3 < 4 ) - sin0|V+> (1 ' 2) , |$ 4 > = cos#|VO (3,4) - sin6»|^_) (1 ' 2) , 

|$ 5 ) = cos#|0 + ) (1 ' 2) + sin6>|0 + >( 3 ' 4 \ |$ 6 ) = cos^|</»_) (1 ' 2) +sin6'| ( /)_) (3 ' 4) , 

|$ 7 ) = cos^|0 + )( 3 ' 4 ) - sin^|0 + )( 1 ' 2 ), |$ 8 ) = cos#|0_) (3 ' 4) - sin^l^ 1 ' 2 ), 

|$ 9 > = cos#|0+> {2 ' 4) - sin 9\<f>+) (1 ' 3) , |$ 10 > = cos^|0 + ) (1 ' 3) + sin#|0 + > (2 ' 4) , 

|$ n ) = cos#|0_) (2 ' 4) - sin#|</>_) (1 ' 3) , |$ 12 ) = cos#|</>_) {1 ' 3) + sinfl|0_) (2 ' 4) , 

|$ 13 > = cosfl|0+> (2 ' 3) - sra0|<£+) (1,4) , |$ 14 > = cos#|0 + > (1 ' 4) + sin #|0 + > (2 ' 3) , 
|$ 15 ) = cos6>|^)( 2 ' 3 ) - sin^l^)^- 4 ), |$ 16 ) = cos#|c/>_) (1 ' 4) +sinfl|0_)( 2 ' 3 ). 



(2.4) 

In this work, we choose 9 = f in the pure states (12.41) . then we obtain the so-called Bell- 
type states which are maximally entangled states. Now, based on momentum distribution we 
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have two types of pure states as the following: 
The first type of pure states 

The first type is with i = add number, this type of pure states are made by Yllj=i ^APiij)® 
\Piij) or Y^ 2 {ijtj}=i MPi'i) ® \P%i3)i which have the same momentum pi for two particles. 
The second type of pure states 

The second type is |$*) with i = even number, this pure states are made by Y^fj=t MPi; ® 
\pj,j) or j=i *) ® bi)*)) which have different momentum for each particle. 
Two types of mixed density matrices 

Consider a two-particle quantum mixed state which is defined as a convex combination of the 
pure states (12.41) . i.e. , 

16 

P = j2q^ i )m, (2.5) 
i=i 

where > 0, Y^] 6 Qi = 1- According to the definition of two types of pure states, we construct 
two types of mixed density matrices as 

16 

Pi = , < % < 1 , = 1, i = odd number, (2.6) 

i=i 

16 

p 2 = , Q < g j < 1 , y^g 7 - = 1, j= even number. (2.7) 

Relativistic quantum states 

In this article, we assume that spin is in the z-direction and momentums are in the yz-plane, 
i.e., p\(2) = (0,pi(2) sin 6*1(2), pi( 2 ) cos 6 l i(2)) and for an observer in another reference frame S 
described by an arbitrary Lorentz boost A, the transformed pure states (see ( lA-iiij) in Appendix 
A ) are given by 

|$*) — > C/(Ai) ® U(A 2 )\&). (2.8) 

For example 

|A$ X ) = U(A 1 )(g)U(A 2 )\^ 1 ) = \Ap u n 1 )\Ap u n 1 ) + \Ap 2 ,n , 2 )\Ap 2 ,n , 2 ) + \Ap 2 ,n 2 )\Ap 2 ,n 2 ) + \Ap 1 ,n' 1 )\A Pu n' 1 ), 
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where 

n m = D{W,p m )\0), n' l{2) = D{W,p m )\l). 



The kets |Api) and \Ap 2 ) are two orthogonal momentum eigen-state after the Lorentz trans- 
formation. Thus, we made two types of mixed density matrices (pi and p 2 ) and after some 
calculation we can see that the mixed state p 2 doesn't change under the Lorentz transforma- 
tion. But, in spacial case, when the momentum and the boost direction are parallel then the 
mixed density matrix p\ which was constructed by the first type of pure states, changes under 
the Lorentz transformation that we will consider in the next section. Let us first obtain the 
feasible region (FR) of the relativistic density matrix p\ as follows: 
Feasible region of the relativistic mixed density matrix pi 

We consider the density matrix of (12. 6p and by imposing PPT conditions with respect to 
each parties, we obtain the FR. For this particular density matrix, the positivity of partial 
transpositions gives linear constraints on the parameters qt (where i is odd number). The PPT 
condition with respect to the first party implies that the eigenvalues of p\ given by 

lefao + g^cos 4 ^ 16(g 12 + g 16 )cos 4 ^ f 

1{2) cos 4 f + cos 4 | ' 3(4) cos 4 | + cos 4 | 



16(g 2 + g5 ) cos 4 16( gl + ge) cos 4 ^ 

A ^ cos 4 1 + cos 4 1 ' 7(8) cos 4 1 + cos 4 1 ' 

16(g 9 - g 13 ) cos 2 | cos 2 | _ 16 (-g 9 + gi 3 ) cos 2 | cos 2 % 

9 " cos 4 % + cos 4 f ' 10 " cos 4 \ + cos 4 f 

16(gii-gi 5 )cos 2 |cos 2 | 16(-g n + g 15 ) cos 2 | cos 2 | 

11 " cos 4 | + cos 4 | ' 12 " cos 4 | + cos 4 f 

16(g 3 -g 5 )cos 2 ! cos 2 | 16(-g3 + g 5 )cos 2 |cos 2 | 

Al3 " cos 4 % + cos 4 1 ' 14 ~ cos 4 % + cos 4 1 ' 

l%i - q 7 ) cos 2 f cos 2 % 16(-gi + q 7 ) cos 2 f cos 2 ^ 

Al5 = cos 4 I + cos 4 I ' Al6 = cos 4 I + cos 4 I ' 
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must be nonnegative, i.e., we have the constraints A» > , for i = 1...16 which give 

99 = 9l3, 9n = 9l5, 93 = 9s, 9i = 97, 

and 

1 

9i + 93 + 9ii +99 = 2' 

where we have used the normalization condition ^\ q i — 1. 

similarly, by imposing the PPT condition to the second party, we obtain 

(~ - 9.) > 0, (2.9) 

It is important to see that the FR of the density matrix pi is independent of the Lorentz 
transformation. 



3 Construction of DOEWs via convex optimization 

We have seen how to manipulate the relativistic density matrix p\ that we want to quantify its 
amount of entanglement, by using DOEW. So the next section deals with the basic definition 
of convex optimization and our scheme to construct DOEW by an exact convex optimization 
method which can be generalized for other density matrices. 

3.1 Entanglement witnesses 

An entanglement witness acting on the Hilbert space 7i = Ti^ £g> Ti^ is a Hermitian operator 
W = ~W\ that satisfies Tr(Wp s ) ^ for any separable state p s in B(7i) (Hilbert space of 
bounded operators), and has at least one negative eigenvalue. If a density matrix p satisfies 
Tr(Wp) < 0, then p is an entangled state and we say that W detects entanglement of the 
density matrix p. The existence of an EW for any entangled state is a direct consequence of 
Hahn-Banach theorem [38] and the fact that the space of separable density operators is convex 
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and closed. Geometrically, EWs can be viewed as hyper planes that separate some entangled 
states from the set of separable states and, hyper plane indicated as a line corresponds to the 
state with Tr[Wp] = 0. 

Definition 1. An EW W is said to be optimal, if for all positive operators V and e > 0, 
the new Hermitian operator 

W = (1 + e)W - eP (3.10) 

is not anymore an EW [39]. Suppose that there is a positive operator V and e > such that 
Wnew = Wopt. — is yet an EW. This means that if Tr(yV op t.p s ) = 0, then Tr{Vp s ) = 0, for 
all separable states p s . Also, one can assume that the positive operator V is a pure projection 
operator, since an arbitrary positive operator can be written as convex combination of pure 
projection operators with positive coefficients. 

Definition 2. When talking about EWs one has to distinguish two different kinds. On 
the one hand, there are the so-called decomposable EWs (DEW), which can be written as 

W = V+Q T l A ) V,Qi>0, (3.11) 

where the operator Q\ is positive semidefmite. It can be easily verified that such witnesses 
cannot detect any bound entangled states. W is non-decomposable EW if it can not be put 
in the form (13. lip (for more details see [10]). One should notice that only non-decomposable 
EWs can detect PPT entangled states. 

3.2 Construction of DOEWs via convex optimization 

At first, we expand an EW for a bipartite system 7i 4 ® 7i 4 as follows 

4 

w = h®h + J2AA®Qj ( 3 - 12 ) 

where J 4 is a 4 x 4 identity matrix, Aij are the parameters whose ranges must be determined 
such that W be DOEW, and QiS ( Q'jS ) are Hermitian operators from the first ( second ) 
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party of the Hilbert space. Via the mapping 

P^TriQ^), P; = Tr(Q' jPs ) (3.13) 
the set of separable states can be viewed as a convex region called FR. 

Let {Q 1 } be bases for the space of Hermitian traceless matrices that operate on 4-dimensional 
Hilbert space that can be written in terms of square matrices E i ., which have the value 1 at 
the position and zeros elsewhere [13], as 

QS = ~^( E l,2 + E 2,l)l Ql = ~^f E U + ^3,1 )> Ql = ^("^V* +" E 4,l)> 

QS = + ^3, 2 )) Ql = ~^- E ^ + ^, 2 )> QS = ^1^3,4 + E 4,s), 

^=■^(^,1-^1,2)) QA = -^( E 4,l - E 1a)i QA = ~^( E 4,2 ~ E 2a)> 
Qa = 7/1^3,1 _ ^1,3)5 QA = 7/1^3,2 - ^2,3); = 7/1^.3 _ ^3,4 )) 

Q 13 = £?!,!, Q M = E 2>2 , Q 15 = E 3 , 3 , Q W = E 4A . (3.14) 

where the subscript S and A indicates symmetric and antisymmetric operators, respectively. 
Notice that for E iy j = where i,j = 1, ...,4, the FR becomes 

4 

E I F mI 2 < 1. (3-15) 
».i=i 

As mentioned above, we will use two steps towards the finding the parameters Aij for the 
density matrix p\ and fully characterize DOEWs based on exact convex optimization method 
for two partite system (see Appendix B). 

The first step 

In the first step, according to the Hermitian traceless basis ( 13.141) . we introduced the maps 
(I3.13P for any separable state p s which map the convex set of separable states to a bounded 
convex region that will be named FR. Then, recalling the definition of an EW, we impose 
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the first condition which is the problem of the minimization of expectation values of witness 
operators W with respect to separable states, i.e., 

rain Tr(Wp s ) > 0, 

where Tr(Wp s ) is the objective function and the inside of the FR will be the inequality 
constraints. So, using the EW of (13.121) after some calculations one arrives at 

4 

Tr(Wp s ) = l+Yl AjPiPj > 0, ( 3 - 16 ) 

which must be satisfied for all Pi and P- belonging to the feasible region. In order to satisfy this 
condition, it is sufficient that the minimum value of Tr(Wp s ) be non-negative. Using standard 
convex optimization, we find this minimum value then we impose the non-negativity condition 
on this minimum. As mentioned before, for the bipartite system the FRs are Y2l=i P? ^ 1 anc ^ 
E]=i Pj 2 < 1 which are in matrix notation as P*P ^ 1 and P n P' ^ 1. Then the Eq. (13.161) can 
be written as a convex optimization problem in the form 

Minimize 1 + P t AP\ 

subject to P*P<1, 
P' l P' ^ 1. 

We now write the primal Lagrangian of this problem, which will be helpful in the following 
development, 

L — 1 + P l AP' - ^A(P'P - 1) - iA'(P /f P' - 1). (3.17) 

where we will assume in the following that A, A' > 0. This problem can be solved using the 
KKT conditions (see Appendix B). From condition VL = , we have 

AP' = AP, P t A=X'P'\ (3.18) 

and from the last KKT condition which we call complementary slackness, we get 

A(P'P-1) = 0, X'iP'P' - 1) = 0, (3.19) 
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which reduce to P t P = 1 and P'*P' = 1. Multiplying AP' = XP by P f from the left side and 
by using P l P = 1, we obtain 

P l AP' = XP f P = A, 
also multiplication P t A = X'P' by P' from the right side and P' P' = 1, yield 

P l AP' = X'P n P' = A', 

then A = A'. Therefore, from ( 13.18P and (I3.19P one can write 

A l AP' = X 2 P', 

and 

AA % P = X 2 P, 

which indicate that A^ for i = 1, 16, are eigenvalues of A 1 A and AA l '. If we choose .4's in 
a way that Af ^ 1 (i.e. — 1 ^ Aj ^ 1) for all z's then the minimum of Lagrangian (13.171) , i.e. 
1 + Al is nonnegative, which leads to the nonnegativity of Tr(Wp s ). 
The second step 

In the second step, for a given p ent , we impose the second condition for an EW, Tr{Wp ent ) < 
0. Now the objective function ( which will be minimized ) is Tr(Wp en t), and the inequality 
constraints come from the first step solution. So, this problem can be written as a convex 
optimization problem of the form 

Minimize Tr(yVp ent ), 

s.t. A 1 A - I d ^ 0, 
We can use associated Lagrange for solving this problem as 

L = Tr(W Pent ) + Tr[(A*A - I d )Z], 
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where Z is a positive symmetric matrix. Then, using the KKT (see Appendix B) conditions 
we get, after some mathematical manipulations, the following formula 

pi + A(Z* + Z) = 0, 

where p\ is a matrix with components 

p^^Tr^.Q^Q 1 .), (3.20) 

A* A = h- 

Obviously, these components are associated with components of witness matrix, A^ v . 
If the matrix Z be a symmetric matrix (i.e. Z l = Z) then p~\ = —2AZ or 

Z=±A% = (3.21) 



where 



A=--piZ~ x . (3.22) 



Now the minimum of the Lagrangian becomes 

l-Tr[y/pWi], (3-23) 

and the negativity of this term is the detection condition for the entanglement of the given 
density matrix. The parameters of matrix A using the Eq. (13.221) are calculated as 

A x ,x = A 2 , 2 = A 3:3 = A AA = - (-= — -i, - 1 , 

2 |6i - 62 | 

A, 2 = A 2 ,i = A 34 = A, 3 = ^( |^_^| - !)» 

-4-6,6 = -4-9,9 = — -4.12,12 = —-4.15,15 = o(t7~ T7 ~~ 

2 |6 5 - 6 6 | 

-4.6,9 = -4.9,6 = —-4-12,15 = — "4.15,12 = ^(jT T7 ~~ -Q) 





-K\ 


h 




\h 


-K\ 




-b 7 



-4-7,7 — -4-8,8 — —-4-13.13 — —-4-14 14 — -( — — - — 1 

2 67 -»8 
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where 

h = 9i + 93 + 95 + 97, ^2 = 99 + qn + 9i3 + 9i5, 

h = qi - q3 - 95 + 97, b 4 = q 9 - q u + gi 3 - g 15 , 

h = qi ~ 93 + 95 - 97, b 6 = q 9 - q u - q 13 + gi 5 , 

&7 = 9i + 93 - 95 - 97, b 8 = q 9 - 9n - 9i3 - 915- 

As an example, let us consider the case b\ > 6 2 , 6 3 > 6 4 , 6 5 > 6 6 and 6 7 > b$ in which 
.Aij's can be written as 

At, j = 1, ^ = 3i 

Ai,j 0, i 7^ j, 



and it is easy to see that 



where 



Fi = g 13 ® g 13 + g 14 ® g 14 + g 15 ® g 15 + g 16 ® g 16 , 
^ 2 = g^ ® g^ + g^ ® g^ - g^ ® g^ - q\ 2 ® q\\ 

Fi = Ql®Ql + Q A s ®Q%-Q\®Q\- Qa ® Q". 
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After some calculations one arrives at 

W = / 4 ®/4-4|$ 1 )(<l> 1 |. (3.24) 

In the next section, we consider the optimality of DEW (13.241) by using Hilbert-Schmidt mea- 
sure and show that the DOEW (13.241) is based on the minimal distance of an entangled state 
from edge of PPT states. 



4 Optimality DEW via Hilbert-Schmidt measure 

According to Ref . [36] . the authors have used the Hilbert-Schmidt distance, which quantifies the 
distance of an entangled state from the set of all separable states and calculated the optimal 
entanglement witnesses explicitly. However, in some situation, specially in high-dimensional 
Hilbert space we are not able to specify the region of separable states. In order to define this 
measure, we recall that the Hilbert-Schmidt norm is defined as 

||4| = V / (A^), (4.25) 

where, (A, B) = Tr(A^B). With help of the norm (14.251) . the Hilbert-Schmidt distance between 
two arbitrary states Pa,Pb can be defined as 

(Ihs(pa,Pb) = \\pa- Pb\\- (4.26) 

By using the Hilbert-Schmidt distance, the so-called Hilbert-Schmidt measure of entanglement 
is defined as 

D(p en t) = min \\p - p ent , || , (4.27) 
pes 

where, S is the set of separable states. In fact, the Hilbert-Schmidt measure is the minimal 
distance of an entangled state p ent . from the set of separable states. 

For an entangled state p en t, the minimum of the Hilbert- Schmidt distance (the Hilbert- 
Schmidt measure) is attained for some state p s since the norm is continuous and the set S is 
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compact. Due to the Bertlmann-Narnhofer-Thirring Theorem [37], there exists an equivalence 
between the Hilbert-Schmidt measure and the concept of optimal entanglement witnesses as 
follows: The Hilbert-Schmidt measure of an entangled state equals the maximal violation of 
the inequality Tr(Wp) > 0, 

D( Pent ) = \\p s ~ P e J\ = -(p ent ,Wopt) = -Tr(p ent W opt ), (4.28) 

where, 

W opt = Ps~P^ -(Ps,Ps- P mt )l ^ (4 2g) 

\\Ps ~ Pent. || 

is an optimal entanglement witness (for more details see Refs. [36, 37J). 

In this work we use this method in order to show the optimality of DEW ( 13.24ft and then 
we discuss the effects of Lorentz transformation on the measure of entanglement of two-particle 
states. One of the important problems is finding nearest separable states to entangled state 
which is not easy for high-dimensional Hilbert space, so we consider the minimal distance of 
an entangled state from edge of PPT states and obtain optimal DEW. 

For an entangled state p ent , the minimum of the Hilbert- Schmidt distance (the Hilbert-Schmidt 
measure) is attained for the edge of PPT states, namely p e d ge - Then we replace p s with p edge 
in (|4T28jl and (T4T29]) as 



V(pent) = Wpedge ~ P ent || = ~ (p ent , V^Dapt) = -Tr( Pent W Dop t) , (4.30) 

and after some mathematical manipulations, we get 



W Dopt = rea9e Pm \. ^'fty = / 4 & J 4 _ 4|$ 1 )($ 1 |, (4.31) 

|| Pedge ~ Pent.\\ 



where p e d ge = Pi according to Eq. (I2.9j) is the edge of PPT states and p ent — pi, so the DEW 
(13T24j) is optimal. 
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5 Entanglement of relativistic two-particle quantum states 

Finally let us consider the amount of entanglement for relativistic pure and mixed states using 
the von Neumann entropy and DOEW, respectively. 

5.1 Entanglement of relativistic pure states 

Let = Y^ij k i=i a ijki\Pi,j) ® \Pki 0> a ijki £ Cbe a general two-particle pure state with nor- 
malization Ylij k i=i \ a ijki\ 2 — 1- F° r this pure state the entanglement of formation E is defined 
as the entropy of either of the two sub-systems, E{\ip)) = —Tr{o\Logi<J\) = —Tr{o<iEog<io<i) 1 
where <i\ (respectively, cr 2 ) is the partial trace of \ip){4>\ over the first (respectively, the second) 
Hilbert space. If Af denotes the matrix with entries given by a^ki, then <j\ can be expressed 
as o\ = NN^ (JV^ means conjugate transpose of M). 

In the previous section we considered two types of pure states which were different in 
combination of momentum eigen-state. If we calculate the von Neumann entropy for the pure 
states (12 .4p under Lorentz transformation with arbitrary boost and momentum direction, we 
will see that these pure states does not change, except for the pure states of the first type 
when the momentum and boost are parallel. 

E(\A&)) — — ^ X k Log 2 Xk, i = odd number, 

k 

Xl(2) 



2(cos 4 (f)+cos 4 (f)) 



cos 4 (|) 

A3{4) - 2(cos 4 (^)+cos 4 (f))' (5 - 32) 



where Xk is eigenvalues of the reduced density matrix (<7i or cr 2 )- After some mathematical 
manipulations, we get 

-1 cos 4 (|) 4 9 2 cos 4 (|) 

£ d* » = cos4( i )+cos 4 ( |){ c08 (2 )to9( 2(co S 4 (f) + co S 4 (f)) )+COS ^MfiMfii' 1 ' 

(5.33) 
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From this expression it can be appreciated that for Q\ = 6 2 when two particles in moving frame 
have the same velocity then the entanglement does not change but for 9\ 7^ 62 the entropy 
decreases with increasing the velocities of the observers. 

5.2 Entanglement of relativistic mixed states using DOEW 

According to Eq. (13.231) and using the complete set of matrices {Q % } for p\ , we get the following 
result: 

Tr[W Pl ] = 1 - Trlv^Tfh] = -(1 - \b x - b 2 \ - \b 3 - 64 1 - \h + h\) 

cos 2 §1 cog 2 62 

-(\h-b 6 \ + \bs + h\ + \b 7 -b 8 \ + \b 7 + b 8 \) 49i 2 f^, (5.34) 

cos 4 ^ + COS ^ 

This result shows that the entanglement between two particles which are described by the 
mixed density matrix pi, under Lorentz transformations is not Lorentz- invariant and decreases 
by increasing the velocity of the observable. 

For detection of entanglement of the relativistic mixed density matrix pi using DOEW (I3.24p 
, the following result is obtained 

cos 2 h. cog 2 <?2 

Tr[W. Pl ] = 1 - 2 Ql - 2q 7 + 4(q 7 - gl ) 2 fe. (5.35) 

cos 4 y + cos 4 y 

Therefore, we have 

Tr[W. Pl ] > Tr[W. Pl } rest . (5.36) 

This result indicates that when momentum and boost are parallel then the entanglement is not 
Lorentz-invariant because Tr[W.pi} res t is more negative than (I5.35p . Therefore, these DOEWs 
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can be used to quantify the amount of entanglement. At the end of this section we show that 
the result (15.351) can be associated with generalized concurrence and coincides with previous 
result for pure state in (I5.33p . 

According to Ref.jH] the entanglement of formation of is given by 

E^Q 1 )) = —n\ 1 Log 2 \i — mXiLogiXi 

where Ai,2 are defined in (15. 32ft and n = m = 2, Ai + A2 = |. Using the expectation value (I5.36P 
for density matrix |$ 1 )(<l )1 | we have 

cos 2 % cos 2 £ r— 
Tr[W. Pl ] = -1 - 4 2 2 = -1 - 8y^\ 2l 

COS 4 + COS ^ 

where q^s are zero except for qi = 1. That is, if we define x — 1 + 8VA1A2, then 



where d is generalized concurrence which has been defined as d = 4a/A]A~2. 

6 Conclusions 

We have obtained DOEWs by using the convex optimization method and shown that it is 
possible to use DOEW operators to quantify the amount of entanglement of bipartite mixed 
states for high-dimensional Hilbert space. This was due to the fact that, one can calculate 
the minimum distance of an entangled state from the edge of bound entangled states, via the 
optimality DEWs. Then, by using it, we have discussed the entanglement of quantum states 
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of two s = 2 fermions under an arbitrary Lorentz transformation. For convenience, instead 
of using the superposition of momenta we used only two momentum eigen states (pi and p 2 ). 
Consequently, in 2D momentum subspace we have shown that for relativistic mixed density 
matrix p\ when momentum and boost are parallel the entanglement is not Lorentz- invariant. 
Our method can be generalized for applying the techniques of this paper to quantifying the 
amount of entanglement of bound entangled states that we can evaluate the minimum distance 
of bound entangle state from the edge of PPT states but in the orientation of inside to separable 
states direction. 



APPENDIX A 

Wigner representation for spin-^: 

It follows [42J that the effect of an arbitrary Lorentz transformation A unitarily implemented 
as U(A) on single-particle states is as follow 



U(A)(\p) ® |<r» = J^^DMW&pMAp) ® (A-i) 

V P a ' 

where 

W(A,p) = L- 1 (Ap)AL{ P ), (A-ii) 



is the Wigner rotation [35] .The Wigner rotation is an element of the spacial rotation group 
SO (3) or subgroup of the homogeneous Lorentz group since it leaves the rest momentum k u 
unchanged : 

W»W = hf. 

This subgroup is called (Wigner's) little group. We will consider two reference frames in this 
work: one is the rest frame S and the other is the moving frame S' in which a particle whose 
four-momentum p in S is seen as boosted with the velocity v. By setting the boost and particle 
moving directions in the rest frame to be v with e as the normal vector in the boost direction 
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and p{(2), respectively, and n = e x pip), the Wigner representation for spin-1/2 is found as 

D*(W(A,p l{2) )) =cos^ + zsin^(a.n), (A-iii) 



where 



fi Pl 7 2) cosh f cosh | + sinh f sinh |(e.j9i7 2 )) . . . 

cos — = — ; (A-iv) 

[| + | cosh a cosh 5 + | sinh a smh.8(e.p{(2))] 



. ^piU . sinhf sinh | (e xpi* (2) ) 

sm — n = — (A-v) 

1 1 + | cosh a cosh 5 + | sinh a sinh 5 (e.pipj))] 



v cos^— ) +(sm— ^n) = 1, 



and 



cosh a = 7 = — , cosh 5 = — , 3 = —. 

^/V^W rn y c 



APPENDIX B 

Convex optimization and Karush-Kuhn- Tucker (KKT) Theorem: 

Many problems from quantum information theory can easily be translated to the language 
of convex optimization problems, of the form 

minimize fo{ x ) 

subject to fi(x) < 0, i = l, ...,m. (A-vi) 

hj(x) = 0, j = l,...,p. 
where x is a vector of decision variables, and the functions f , /j and hj , respectively, are the 

cost, inequality and equality constrains which satisfy inequality fi(ax + By) ^ afi(x) +3fi(y), 
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for all x,y G R and all a, (3 G R with a + (3 = 1, a ^ 0, /9^0 because of convexity and 
the equality constraint functions hi(x) = must be affine (A set C G R n is affine if the line 
through any two distinct points in C lies in C). 

There are a number of important necessary conditions that hold for problems with zero dual- 
ity gap. These Karush-Kuhn- Tucker conditions turn out to be sufficient for convex optimiza- 
tion problems. Assuming that functions and hj are differentiable and that strong duality 
holds, there exists vectors ( e i? fc , and y G R m , such that the gradient of dual Lagrangian 
L(x*, (*, y*) = f(x*) + J2i Qhi(x*) + J2i U*i9i{ x *) over x vanishes at x* 

hi(x*) = (primal feasible) 
9i(x*) < (primal f easible) 
y* > (dual f easible) 

y*gi(x*) = o 



Then x* and ((*,y*) are primal and dual optimal, with zero duality gap. In summary, for 
any convex optimization problem with differentiable objective and constraint functions, any 
points that satisfy the KKT conditions are primal and dual optimal, and have zero duality 
gap. Necessary KKT conditions satisfied by any primal and dual optimal pair and for convex 
problems, KKT conditions are also sufficient (see Ref. p[5j for more details) . 
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